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ABSTRACT 


An  algorithm  for  designing  non-recursive  frequency- 
selective  digital  filters  having  equal  stop-band  and 
pass-band  ripples  of  specified  magnitude  is  described. 
There  is  no  known  upper  limit  on  the  order  of  the  filter 
that  can  be  designed  by  this  algorithm  and  filters  with 
orders  in  the  hundreds  already  have  been  designed. 
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A  NEW  TECHNIQUE  FOR  THE 
DESIGN  OF  NON-RECURSIVE  DIGITAL  FILTERS 


* 


The  purpose  of  this  note  is  to  give  a  preliminary  description  of  some  new 
results  in  the  design  of  non-recursive,  digital  filters.  The  specific  filter  types 
under  consideration  are  low-pass,  band-pass,  band-reject,  etc. ,  although  there 
are  reasons  to  believe  that  the  technique  to  be  described  can  be  extended  to  more 
general  types  of  filters. 


An  excellent  survey  of  the  state-of-the-art  of  non-recursive  filter  design  is 

given  in  a  recent  paper  by  Rabiner  .  Another  useful  reference  in  this  connection 

2 

is  Rabiner,  Gold  and  McGonegal  .  The  present  discussion  is  concerned  with  a 

much  improved  method  of  designing  linear-phase,  equi-ripple  filters  of  the  type 
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proposed  by  Herrmann  and  Schuessler  ’  .  These  filters  are  designed  to  have  a 

specified  number  of  ripples  in  the  pass -band  (N  )  and  in  the  stop-band  (N  ).  In 

P  s 

addition,  the  magnitude  of  the  ripple  is  specified  in  the  pass-band  (1  ±  6  and  in 
the  stop-band  (±62)*  The  equi-ripple  filter  meeting  these  requirements  is  one 
whose  frequency  response  has  N  +  N  extrema  that  achieve  the  maximum  allowable 

P  s  j  2nf T 

ripple  values.  An  example  is  shown  in  Fig.  1.  In  this  figure,  H  (eJ  )  is  the 
frequency  response  of  the  filter  at  frequency  f  and  T  is  the  time-domain  sampling 
interval. 


Since  the  frequency  response  of  a  non-recursive  digital  filter  is  a  trigono¬ 
metric  polynomial,  Herrmann  and  Schuessler  attempted  to  solve  the  equi-ripple 

filter  design  problem  by  writing  a  set  of  2  (N  +  N  -  1)  non-linear  equations  that 

th  ^  ^ 

force  an  (N^  +  Ng  -  1)  order  polynomial  to  have  extrema  that  achieve  the  maximum 
allowable  ripple  values.  These  equations  are  then  solved  by  means  of  a  non-linear 
programming  technique  which,  unfortunately,  fails  to  converge  when  the  order  of 
the  polynomial  is  larger  than  about  twenty.  This  is  the  main  drawback  to  their 
method.  The  design  algorithm  that  now  will  be  described  has  no  known  limitation 
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Fig.  1.  An  example  of  an  equi-ripple  filter. 
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on  order  and  filters  having  orders  in  the  hundreds  already  have  been  designed  with 
it. 


The  proposed  design  algorithm  is  basically  an  iterative  technique  for  pro¬ 
ducing  a  polynomial  that  has  extrema  of  preassigned  value.  It  begins  by  making 
an  initial  estimate  of  the  frequencies  at  which  these  extrema  will  occur  and  then 
uses  the  Lagrange  interpolation  formula^  to  obtain  a  polynomial  that  goes  through 
the  maximum  allowable  ripple  values  at  these  frequencies.  The  details  of  this 
initial  guess  are  not  important;  the  author  has  been  using  a  uniformly  spaced  set 
of  frequencies  ranging  from  0  to  0.5  and  requiring  the  polynomial  to  alternatingly 
go  through  1  ±  6  j  in  the  pass-band  and  ±  62  in  the  stop-band.  This  initial  set  of 
frequencies  along  with  the  associated  Lagrange  interpolation  polynomial  are  sketched 
in  Fig.  2.  Note  that  the  polynomial  associated  with  the  initial  guess  does  not  have 
extrema  that  achieve  the  maximum  allowable  ripples  but,  rather,  extrema  that 
exceed  these  values.  The  next  stage  of  the  algorithm  is  to  locate  the  frequencies 
at  which  the  extrema  of  the  first  Lagrange  interpolation  polynomial  occur.  These 
frequencies  are  taken  to  be  a  second,  hopefully  improved,  guess  as  to  the  frequen¬ 
cies  at  which  the  extrema  of  the  filter  response  will  achieve  the  desired  ripple 
values.  This  second  set  of  frequencies  is  indicated  in  Fig.  2.  The  algorithm 
now  "closes  the  loop"  by  using  these  new  frequencies  to  construct  a  Lagrange 
interpolation  polynomial  that  achieves  the  desired  ripple  values  at  these  frequen¬ 
cies.  The  extrema  of  this  new  polynomial  are  then  located  and  used  to  start  the 
next  cycle  of  the  algorithm.  The  algorithm  is  reminiscent  of,  but  different  from, 
the  Remes  exchange  algorithm  used  in  the  theory  of  Tchebycheff  approximation. 

The  formal,  mathematical  description  of  the  design  algorithm  can  be  written 

til 

as  follows.  Denote  the  frequencies  obtained  at  the  i  stage  of  the  algorithm  by 
f  jj),  k  =  1,  2. . .  N  +  Ng,  f^  =  0,  f^  +  =  1/2  T.  The  frequency  response  of 

*  n  G 
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Fig.  2.  The  design  algorithm  for  N  =  N  =3. 
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the  filter  at  the  ith  stage  of  the  algorithm  is  then  given  by 


G(l)  (f)  = 


N  +  N 
P  s 

Z 

k  =  1 


N  +  N 
P  s 

Z 

k  =  1 


A«  Y, 
k  k 

(X  -  X^) 


l(0 


(X  - 


where 


X  =  cos  2nfT,  x£*  =  cos  2nf^  T,  G^  (f)  =  H^l)  (e-*2"^) 


,0)  - 


N  +  N 
P  s 

n 

e  =  1 
e/k 


,-l 


(X^  -  X(i)) 
k  e 


and  Y  =  1  ±  6  .  for  subscripts  referring  to  the  pass-band  and  Y  =  ±  6  for  subscripts 

rC  1  K  Z  ^ 

referring  to  the  stop-band.  The  frequencies  associated  with  the  (i  +  1)  stage  of  the 
algorithm  are  now  defined  to  be  the  +  Ng  solutions  to  the  equations 


d  G(l)  (fk*l+1)) 
df 


=  0  ,  k=l,. 


N  +  N 
P  s 


defining  the  N  +  N  extrema  of  G^(f).  The  (i+l)S*  frequency  response  G^+^(f) 

p  S  /.\ 

can  now  be  obtained  by  replacing  X  ^  by  X  k  in  the  above  equation.  This  process 
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is  then  repeated  until  the  desired  degree  of  convergence  has  been  obtained. 

It  should  be  obvious  that,  if  the  procedure  just  described  converges  it  must 

yield  the  desired  equi-ripple  filter.  The  algorithm  has  been  programmed  by  the 

author  on  a  Hewlett-Packard  9100A  desk  calculator  and  on  the  I.B.  M.  360  computer 

by  J.  Siegal,  an  M.  I.T.  graduate  student.  On  both  machines  the  algorithm  has  been 

demonstrated  to  converge  to  within  a  preset  accuracy  limit  in  a  surprisingly  small 

st 

number  of  iterations.  An  example  of  a  41  order  band -pass  filter  designed  by  the 
Hewlett-Packard  program  is  shown  in  Fig.  3  and  an  example  of  a  251  order  low- 
pass  filter  designed  by  the  I.B.  M.  360  program  is  shown  in  Fig.  4.  The  latter 
design  required  12  iterations. 

There  are  a  number  of  interesting  mathematical  conjectures  that  can  be  made 
about  equi-ripple  filters  and  the  above  described  algorithm  for  designing  them. 

The  author  in  conjunction  with  Prof.  A.  Oppenheim  of  MIT  believe  they  can  prove 
these  conjectures;  however,  full  details  are  not  available  at  this  time.  They  will 
be  published  as  a  paper  elsewhere  in  the  literature  in  the  near  future.  Because  of 
the  pertinence  of  these  conjectures  to  the  above  algorithm,  it  seems  appropriate 
to  state  the  more  important  of  them  here.  They  are  as  follows: 

1.  The  equi-ripple  filter  is  optimum  in  the  sense  that  no  other  non¬ 
recursive  filter  having  the  same  number  of  pass-band  and  stop-band  ripples  and 
whose  ripples  are  no  greater  than  1  ±  6  ^  and  ±  ^  *n  the  pass -band  and  stop-band 
respectively  can  have  a  smaller  transit  ion -band.  (The  transition -band  is  defined 
as  that  interval  of  frequencies  for  which  62<H(f)<l-6j.) 

2.  The  algorithm  described  above  always  converges  to  that  unique  equi- 
ripple  filter  meeting  the  pass-band  and  stop-band  specifications. 
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Fig.  3.  A  41st  order  band -pass  filter  designed  for  Np  =  9, 
N  =  12,  6,  =  .005,  60  =  10 "5. 

Sl  2 
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Fig.  4.  A  251st  order  low-pass  filter  designed  for  N  =  32,  N  =  94,  6.  =  .01,  6„  =  4  X  10  \ 

p  S  X  z 
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